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In non-minimal Higgs sectors, coupling constants of the discovered Higgs boson can deviate from 
the predictions in the Standard Model by effects of additional scalar bosons. The pattern of the de¬ 
viations in various Higgs boson couplings largely depends on the structure of extended Higgs sectors. 
Therefore, we may be able to determine the true Higgs sector by hngerprinting the predictions on 
the Higgs boson couplings in each model with the future precision data. As the expected precision 
is extremely high, it is essentially important to evaluate the theory predictions as accurately as 
possible with radiative corrections in each model. We calculate a full set of one-loop corrections 
to the Higgs boson couplings in the two Higgs doublet models with a softly-broken Z 2 symmetry, 
where there are four types of Yukawa interactions. We discuss how to distinguish four types of the 
model by evaluating the pattern of the deviations in the Yukawa coupling constants and also the 
couplings with gauge bosons. In addition, we demonstrate how inner parameters in each type of 
the model can be extracted by the future precision measurements of these couplings at the high 
luminosity LHC and the International Linear Collider. 


I. INTRODUCTION 

The LHC Run-I has led the Standard Model (SM) particle content to be completed by discovering the Higgs 
boson (h) with the 125 GeV mass ii- The data tell us that properties of the Higgs boson are consistent with 
those of the Higgs boson in the SM within the current uncertainties. 

The discovery of h does not directly mean that the true Higgs sector is the SM Higgs sector. In the SM, 
the Higgs sector is assumed to be the minimum form to spontaneously break the SU{2)l x U{1)y gauge 
symmetry into t/(l)EM- There are possibilities for extended structures of the Higgs sector with some number 
and representation of additional Higgs fields; e.g., models with the doublet field and additional singlet fields, 
doublet fields and/or triplet fields. These extended Higgs sectors, like the SM Higgs sector, can also explain the 
current LHC data for h in some portions of their parameter regions. Therefore it is important to determine the 
true Higgs sector by future experiments. 

Many of new physics models beyond the SM have characteristic extended Higgs sectors. For example, by 
introducing additional Higgs fields, new physics phenomena such as neutrino mass, baryon asymmetry of the 
universe and dark matter can be explained. There is a strong relation between the structure of the Higgs sector 
and the new physics theory. When the true Higgs sector is determined by future experiments, we can obtain 
an important hint for new physics. Therefore, the study on the determination of the Higgs sector is one of the 
most important studies of new physics. 

It has been known that indirect tests for new physics can be performed by utilizing precision measurements 
of electroweak observables such as the electroweak oblique parameters at the LEP/SLC experiments. By the 
discovery of the Higgs boson and the measurement of their properties, the Higgs boson couplings have been added 
to the list of such useful observables. It is important to study deviations of the Higgs boson couplings from the 
SM predictions in various new physics models, because measurement accuracy of these coupling constants will 
be improved at future experiments such as the High-Luminosity LHC (HL-LHC) and future electron-positron 
collider experiments. At the International Linear Collider (ILC), most of the Higgs couplings are expected to 
be measured with typically 1 % or better accuracy M- In order to use these future precision measurements, 
we should calculate various Higgs couplings as precisely as possible with radiative corrections. Furthermore, by 
evaluating the pattern of deviations in the Higgs couplings in each extended Higgs sector, we may be able to 
distinguish extended Higgs sector. 

There have already been many studies on radiative corrections to the h couplings in extended Higgs sectors. 
In Refs. [M3, onc-loop corrections to the h couplings have been studied in two Higgs doublet models (THDMs). 
Radiative corrections to the h couplings in the Higgs triplet model have been calculated in Refs. fTTI |. 

In this talk, we investigate radiative corrections to the h couplings in THDMs with a softly-broken Z 2 
symmetry as an example of extended Higgs sectors. Under the Z 2 symmetry, four types of Yukawa interactions 
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TABLE L Charge assignment of the Z 2 symmetry given in Ref. [l^ and mixing factors in the Yukawa couplings. 


can appear depending on the assignment of the Z 2 charge to particles. We calculate renormalized couplings 
hVV, h^ and hhh {V = W and Z, and f = c, t, b and r) at the one-loop level in the on-shell renormalization 
scheme Q . We discuss how to distinguish four types of the THDM by evaluating the pattern of the deviations in 
the Yukawa coupling constants and also the couplings with gauge bosons. We further numerically demonstrate 
how inner parameters of the models, such as masses of extra Higgs bosons and mixing parameters, can be 
extracted when scale factors are experimentally determined with expected uncertainties at the HL-LHC and 
the ILC. We take into account current experimental bounds and theoretical bounds from perturbative unitarity 
and vacuum stability in our analysis. 


II. TWO HIGGS DOUBLET MODELS 


In the THDMs, there are two iso-spin doublet Higgs fields <i>i and ^ 2 - ‘f’l and $2 rece ive vacu um expectation 
values (VEVs) vi and V 2 , respectively, where they satisfy relations u(~ 246GeV) = -I- 

and tan/3 = V 2 fvi. After the electroweak symmetry breaking, five physical mass eigenstates (i.e., charged 
Higgs bosons a CP-odd Higgs boson A and two CP-even Higgs bosons h, H) and three unphysical Nambu- 
Coldstone bosons and appear. 

We impose a softly-broken Z 2 symmetry in order to forbid flavor changing neutral currents at the tree level. 
There can be four different types of Yukawa interaction, depending on the assignment for the charge of the Z 2 
symmetry. We define the Z 2 charges for all the particles as shown in Table I [Tj, 0] • 

We assume CP invariance in our present analysis. We can write the Higgs potential as 


U = + ml\^2? - ml{^\^2 + + yI$ 2|^ + A3|$in$2| 

+ A4|$I$2p + Y 


( 1 ) 


where all parameters (namely - m|, Ai - A 5 ) are real. m| is the parameter which indicates the soft-breaking 
scale for the Z 2 symmetry. These parameters are replaced by the physical parameters; namely, the masses 
of H^,A,H and h, the two mixing angles a and /3 which correspond to those among CP-even Higgs fields 
and charged (and CP-odd) Higgs fields, respectively, the VEV v and the remaining parameter M^, where 

[|]. 

Sin p cos p « 


THDM 

We here explain the scale factors for SM-like Higgs boson couplings, kx = , where X is any field 

9hxx 

interacting with h. At the tree level, the scale factors of gauge interactions (kv) correspond to sin(/3 — a). 
We define that the limit in sin(/3 — a) approaching to unity is the SM-like limit [2. The scale factors of hff 
coupling constants can be expressed by k/ = sin(/3 — a) + ff cos(/3 — a), where ff are given in Table ID The 
formulae of Kf, and Kt are different in each type of the THDMs. In Ref. [l^, the possibility of fingerprinting 
the THDMs with future data at the HL-LHC and the ILC has been studied in detail by using the tree level 
analysis. The full analysis with radiative corrections have been studied in Refs. [H, 0. 


III. ONE-LOOP CORRECTED HIGGS BOSON COUPLINGS 

In this section, we give the results of our simulations analyses for extraction of inner parameters by hngerprints 
of scaling factors. By the on-shell renormalization scheme, we calculate several scale factors at the one-loop 
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FIG. 1: Scatter plots for a bench mark set {Any, Anb, Ak,t) = (—2%,+5%,+5%). Blue (red) points satisfy the bench 
mark set within the 1-cr uncertainty at the HL-LHC (ILC500) given in Ref. 


level which are defined following as, 


_ f/ixjc[P i,P2j9^]thdm 

Kx = —^- 

^hxx[pl,pl,q^]sM 


( 2 ) 


where 9^]sm(thdm) are renormalized coupling constants in the SM (THDMs). Details of the 

renormalization calculations are described in Ref.[l|. 

We here give approximate formulae of one-loop corrected Akx dehned as Akx = kx — 1 because it is useful 
for understanding our numerical results. We here show approximate formulae of Aky, Akr and Akc, 


Aky — — 
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( 3 ) 


where x = a — (3 + ^ and c$ = 2(1) for d) = H^{A and H). We can see that there are loop contributions 
composed of m|,/t!^(l — M^)^ in Akx- The contributions come from the wave function renormalization for the 
h field. The magnitude of the loop correction strongly depends on the balance of M and v. If \M\ is comparable 
to V, the contribution gives a quadratic power-like dependence as m\. On the other hand, if |M| is much larger 
than V, the loop effect in Aky reduces as l/m| according to the decoupling theorem Akf have the same 
loop contributions because the wave function renormalization appears in all renormalized Higgs couplings. For 
Akb and Akt, because there are non-negligible loop corrections which come from top quark loop diagrams, their 
approximate formulae are not simple as that of Akr [H- 

We suppose that Akv, Anb and Ak,- will be measured at the HL-LHC and the ILC500 with uncertain¬ 
ties expected in Refs. m El- According to Ref. measurement uncertainties of Akx are given as 
[a{Ky),a{K,b),(j{Kr)] = [2%,4%, 2%] for the HL-LHC and [a{Ky),a{Kb),<r{nT)] = [0.4%, 0.9%, 1.9%] for the 
ILC500. We consider a benchmark set for central values of {Any, Ant, AKr) = {—2%,+5%,+5%) as an ex¬ 
ample. We here assume measured Akc to be a negative value, so that the bench mark set values indicate 
that the Higgs sector seems that of the Type-H THDM. We study parameter regions in which values of kx 
are predicted around the central values within the l-cr error described above by scanning inner parameters 
X, tan 13, m$(= m^± = mn = itia) and in the Type-H THDM. In these analyses, we take into account con¬ 
straints from perturbative unitarity and vacuum stability. Considering current experimental bounds in addition 
to the theoretical bounds, we take tan ,5 and as tan ,5 > 1 and to<i> > 300GeV, respectively. 

We demonstrate to extract inner parameters x, tan/3, m$, and C by using future precision data at the 
LH-LHC and the ILC500, where we define such as C = 1 ~ and rh<i, = ( parameter indicates 

the magnitude of non-decouplingness of loop effects of extra Higgs bosons. When f is unity, is equal to 
the mass of extra Higgs bosons m$. If C —>■ 0, namely if the magnitude of closes to that of m|,, the loop 
corrections of extra Higgs bosons become zero. 

Let us start to explain the left panel in Fig. [1] which shows the relation between x and tan /3 in the favored 
regions from assumed future data. Blue points (red points) indicate arrowed regions at the HL-LHC (ILC500). 
We can understand the behaviors by using the approximate formula of Aky — Akr eq. ©. Because values 
are assigned in Aky and Akr, x is inversely proportional to tan/3. As you can see, points follow the inversely 
proportional relation. In the case using uncertainties at the LH-LHC, the arrowed region of Any is —4% < 
AKy iS 0%, which includes the SM-like limit. Because the value of AKy includes the SM-like limit, the limit 
a; —>■ 0 is still allowed. It means that values of mixing angles cannot be determined with precisely at the HL-LHC. 
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On the other hand, at the ILC500, Akv can be constrained to the narrow region such as —2.4% < < —1.6%. 

We can extract values of the mixing angles more precisely as —0.07 ^ x ^ —0.05 and 1 < tan/3 < 1.8. 

The middle panel of Fig. [T]shows the relation between x and m$. At the ILC500, regions with x ~ —0.06 are 
allowed as described in the previous paragraph. m<i> is expected to be non-zero as 400GeV < < 550GeV. 

That is based on the following reasons, x and tan,5 are fixed to be a; ~ —0.05 and tan/3 ~ 1, respectively, in 
order to satisfy Akt-({,) ~ —5%. Because it is impossible to be the non-zero value of Aky by only the mixing 
effects, non-zero loop corrections are required. On the other hand, at the HL-LHC, rh$ can be still zero, because 
the region Aky —>■ 0 is allowed. 

Finally, we discuss the behavior of the third panel in Fig. [T] , which shows the m$-C plane. We can find 
the upper bound for the mass of extra Higgs bosons; i.e., < SOOGeV at the ILG500. There are also upper 

and lower bounds for f for each value of m$. We can extract non-zero values of C as extraction of fhii, by 
using measurement uncertainties at the ILC500; i.e., 0.7 ^ C ~ 1-6- Therefore, the non-decouplings effect of 
additional scalar bosons can be extracted. Such information is useful to narrow down models of new physics. 
On the other hand, at the HL-LHG, regions with C = 0 are allowed for all values of m,s,. 


IV. CONCLUSION 

We investigated radiative corrections to the discovered Higgs boson couplings in all types of THDMs with 
the softly-broken Z 2 symmetry at the one-loop level. We demonstrated how precisely we can extract inner 
parameters of extended Higgs sectors by comparing our theoretical calculations of the Higgs boson couplings 
with future precision measurements at the HL-LHC and the ILC500. In our analysis, we considered theoretical 
constraints such as perturbative unitarity and vacuum stability as well as current experimental data. We 
found that mixing parameters such as x and tan /3 can be determined more precisely by using measurement 
uncertainties at the ILC. Furthermore, there are possibilities to obtain the upper bound for the mass of extra 
Higgs bosons without their direct discoveries and also to get information of the decoupling property. The 
determination of the Higgs sector is one of the most important studies of new physics. In order to determine the 
structure of the Higgs sector by fingerprinting the Higgs boson couplings, the comprehensive study of radiative 
corrections to the Higgs boson couplings is an important task. 
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